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Using the expansion of the baryon wave function in a series of products of single 
quark bispinors (Dirac orbitals) , the nonsinglet axial and tensor charges of a nucleon 
I are calculated. The leading term yields gA = 1-27 in good agreement with exper- 

iment. Calculation is essentially parameter- free and depends only on the strong 

■ coupling constant value q^. The importance of lower Dirac bispinor component, 
ijT) ', yielding 18% to the wave function normalization is stressed. As a check, the baryon 
^ I decuplet masses in the formalism of this model are also computed using standard 

■ values of the string tension a and the strange quark mass m^; the results being in 
, a good agreement with experiment. 
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Axial and tensor charges of a nucleon are important to characterize the basic structure 
of the nucleon as composed of strongly coupled quarks PP^Il]. Whereas many baryon 
charateristics can be reasonably obtained in relativized quark models where all relativis- 
tic effects are treated via Salpeter equation and spin corrections, the axial and tensor 
charges are sensitive to the Dirac structure of quark wave functions, and in particular, as 
will be shown below, to the negative energy components. The most systematic treatment 
of relativistic baryon system can be done in the three-body Bethe-Salpeter formalism, 
however in this approach a rigorous derivation ends up in the system of more than 20 
integral equations, and therefore, a drastic simplification is needed to make realistic 
calculations. Recently this kind of approach was developed in [Hj, where quark-quark in- 
teractions have been properly parametrised. However the relativistic objects like negative 
energy component admixture rj (see below) are very sensitive to the form of interactions, 
and in jHj ?7 has come out rather small. In our paper we choose another and much simpler 
approach, which allows to treat all Dirac components properly including lower ones, and 
to work out the resulting rj (and axial and tensor charges) in a transparent way. 

First of all, let us recall some formulas of the spin relativistic theory. Consider a 
polarized free fermion with spin |. The fermion polarization could be described by a 
spacelike pseudovector a^, which reduces in the fermion rest frame to the unit 3- vector 
along which the fermion spin projection is equal to +1/2. The pseudovector obviously 
meets conditions: a ■ a = —1, a ■ p = 0, where p is the the fermion 4- momentum. The 
corresponding Michel-Wighimann fermion density matrix is 

(1) 

where m is the fermion mass. Note, that here and below we propose the plane wave 
amplitude being normalized to unity: uu = 1. Using the density matrix the pseu- 
dovector could be also defined in an arbitrary frame as the averaged value of the 
Pauli-Lubansky operator: 

(W^^) = la^ ly^ = --l-eA^-P-X^P^, e°'^' = +1. (2) 

where J^^ is the angular momentum operator. 

Axial and tensor current operators are defined as follows: 

= V'TmT'V', jj. = (3) 
where cr^ = ^[7/i,7i/] (both currents are hermitean). It it easy to check^ that for a plane 

(j^) = w7m7'w = (ij.) = ^a^.7'M = - «V) (4) 

Now, consider a polarized nucleon, e.g. proton, with momentum p and polarization 
a, and define two quantities — the axial charge qa and the tensor charge qt in terms of 
average values of quark axial and tensor currents similarly to (^, as follows: 

(p, a u'jf.'y^u - d-ff,-f^d p, = -gAa^i, (5) 



wave 



^The simplest way to derive these relations is to work in the fermion rest frame and to choose the Oz 
axis direction being parallel to the space part of the pseudovector a 
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p,a 



p^a) = —(a^'p'' -ayn (6) 
m 



where m is the proton mass. 



In the hmit of exact isotopic symmetry (which we assume further), the isotopic struc- 
ture of the axial current is 2/3, and its matrix element must coincide with that one of 
the current from the same isotopic multiplet with structure /+ taken between proton 
and neutron. The latter transition arises in neutron /5-decay and we conclude that the 
constant qa must coincide with the neutron decay axial constant Ga, which is very well 
known experimentally^ Ga = 1.27. For the tensor charge, unfortunately, we have no well 
estabhshed experimental value, although it could be in principal measured as the first 
momentum of the proton transversity distribution function hi{x) [Zj. 

As for the theoretical calculations of the charges, it is well known that the nonrela- 
tivistic quark model predicts qa = ^, while for massless relativistic quarks in the MIT 
bag model, one obtains much smaller values qa = 1.09, both in disagreement with exper- 
iment. As for the tensor charge gx, various theoretical estimations give disparate results 
ranging from 0.89 to 1.45, see [Hj for Refs. and discussions. All this certainly stimulates 
further study of the question. 

The approach proposed in this paper could be applied not only for analysis of nucleon 
properties, but also for the calculation of constants of hyperon semileptonic decays, which 
values are quite important for the accurate determination of the CKM matrix elements 
(see, e.g., 0). We plan to develop a systematic analysis of baryon properties in the 
framework of Dirac orbital expansion in future papers; the first part of which is reported 
below. 

To construct the baryon wave function, one starts with the Hamiltonian ^0] obtained 
in the instantaneous approximation from the general Bethe-Salpeter equation: 

3 

^$(ri,r2,r3) =E$, H = Y,H, + AH (7) 

i=l 

with 

Hi = p(i)Q:(i) + /5(i)(mi + M{ri)) (8) 
where M(rj) in the limit of vanishing gluon correlation length is 

M = a\ri-rx\ (9) 

a is the string tension, and rx is the string-junction coordinate. Here AH contains 
perturbative gluon exchanges. We expand the baryon wave function in a series of products 
of quark eigenfunctions ipn^ = (^{i)) , namely [TT] 

3 

vl>(ri,r2,r3) = J^H^S {r.)Cn,n,n, (10) 

{rii} i=l 



In what follows we shall consider the leading valence approximation for the nucleon 
keeping only the first term in (fTUj) . \E' — \l/o, which contains the ground state S'-wave 

^Our choice of the Ga sign corresponds to 
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Dirac orbitals \ux) and \dx) for u and d quarks with spins up and down. So, for the 
proton polarized along the axis Oz, one has: 

\P T) = [-2(|m ^uUi)+ perm.) + {\u ^ i d ^) + perm.)] (11) 

This expression has the same form as in the standard SU (4) or SU{6) model |3] except 
for the bispinor contents of \ux) and I^a)- In (fTTj) ux and dx orbitals depends on space 
coordinates via the only function x\ in accordance with the isotopic symmetry principle. 
Note, the proton state (|TT| with S-wave orbitals evidently has the average 3-momentum 
being equal to zero: 

(P) = J2(P^) = (12) 

where pj is the 3-momentum of the i^^ quark. 
Insertion (|TT|l into © and (jHl) yields^: 

9a = (p^ n+SgU - rf+Sgc/ p t) = +^(XtI^3|Xt) - ^(Xil^slXi) = +^(XtI^3|Xt) (13) 



(14) 

where we use obvious equalities: 



(XtI^sIxt) = -(Xil^slXi), {xi\P^3\xd = -{xi\P^3\xi) (15) 

To proceed further, we choose the standard representation of the bispinor wave func- 
tion 

°^ 13= (I °] (16) 



aay ' \0 -1 

, , 1/G(r)fi.,o,i(^,0)\ 
"^^^''''^^Kf(.)44(M)J ^^'^ 



The following normalization conditions hold: 

oo 

J nj,i,^,njirndo = 5jf5ii'5^rn', j [G'^(r) + F^(r)] dr = 1 



(19) 



To take into account perturbative gluon exchanges we represent AH effectively as 
one-particle operators 



^see for similar calculations in bag models 
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Table 1: gA, Qt, and r] for various theoretical prescriptions in comparison with experi- 
mental data 





Exp. 


NRQM 


C = 


C = 0.3 


9A 


1.27 


1.67 


1.36 


1.27 


9t 




1.67 


1.51 


1.47 


V 







0.14 


0.18 



with an effective coupling constant ( being fixed by realistic ct^ values (see below), and 
the equations for G{r), F{r) acquire the form [T!^ 



G' - -G - ie + m + ar + -] F = 0, 
r \ r 
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F' + -F+\e-m-ar + -\ G = Q 



(21) 



where m is the light quark mass and e is its energy. 

Now, using equations ()17pi8|l . one can calculate explicitly the matrix elements ap- 
peared in (|13n4jl : 



(xtI^sIxt) 



^G^ir) 



r drdo 



F^(r) 



(fiioilfTsl^ioi)"^ ^(^i 1 i|o"3|^i 1 i) 

\ 2'"'2 2'"'2 r'S 2' '2 2' '2 



dr ( G\r) - -F^{r) 



[22) 



(xtI/^^sIxt) 



r drdo 



GHr) 



1 1 i|o-3|rii 1 1] 

2 ' ' 2 2 ' ' 2 



dr ( G\r) + -F\r] 



(23) 



Finally qa and qt can be written as^ 



9a 



9t 



+1 - 1" 



F^{r)dr 



(24) 



Note, that in the nonrelativistic limit 77 = and qa = 9t = 

The light quark mass m is very small and can be omitted from (j?H) with good accuracy. 
So, the only dimensionful parameter a remains in the model, but both qa and qt constants 
are dimensionless and, correspondingly, do not depend on a. Thus, ( proves to be the 
only essential model parameter. We have computed the values of rj, qa, and qt for two 

*It is clear from the aforesaid that neutron charges are simply opposite in sign to the proton ones, so 
we do not need to consider them separately 
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different values of (: C = ^iiid C = 0.3. The results are given in Tabled there also given 
for comparison the values obtained in the nonrelativistic quark model. 

One can see that the resulting is in the correct ballpark for ( G [0, 0.3]. The value 
( = 0.3 corresponds to the reasonable effective value of in the qq potential, namely 

from (j2 = (Y1 and (rij) ^ V3{ri), one has 

(a.)eff = ^ 0.39 

and this value of was checked in the actual calculation of the baryon masses in the 
Effective Hamiltonian approach ^3]. It is rewarding that the resulting = 1.27 is in 
close agreement with experiment. 

Concerning the tensor charge gx, one can compare these results with the lattice 
data pni, where both gA and gr are close to each other and are in the interval 
1.12 < gA, gr < 1.18 for > 0.5 GeV. 

There are two possible unaccounted effects which can influence our results. First, 
the contribution of other terms in ()10|1 - excited Dirac orbitals. The corresponding 
multichannel calculations done in ^1] for magnetic moments, result in decreasing of 
the modulus of magnetic moments of proton and neutron by some 10 — 15% when one 
accounts for four Dirac orbitals for each quarks, and one can expect the same type of 
corrections for gA- Second, the contribution of chiral degrees of freedom, i.e. of the vr, 
7], K exchanges. Again, for nucleon magnetic moments these corrections are typically of 
the order of 10% [TT], and we expect this to be an upper limit for gA, since magnetic 
moments are much more sensitive to the contribution of the lowest Dirac components, 
than gA and gT-, where these contributions enter quadratically and not linearly. 

These corrections are not taken into account above, which is planned for a subsequent 
work, where also hyperon semileptonic decays are considered [T6] . 

One should note, that relativistic approach to the baryon wave function, based on 
the light-cone formalism jl7] was shown to improve the qualitative agreement of gA with 
experiment, however quantitatively still far from experiment. 

What is also possible to calculate in this approach are the baryon decuplet masses. We 
restrict ourselves here only to the ground state, so our leading S'-wave approximation 
looks rather reasonable. As the spin and the isospin of any pair of light quarks in a 
decuplet baryon is equal to unity, one may expect the pion exchange effects are suppressed 
for decuplet masses. Indeed, calculations in jTTj show that pion exchanges yield a negative 
shift to mA by about 30 MeV. So, the baryon mass is completely determined by the 
eigenvalues of the system of equations ^2,\\ . The only modification is due to the existence 
of a valence s-quark in a baryon - in fact, we should calculate independently the energy 
£ and orbital x fo^^ the light quark q and for the strange quark s; those are obtained from 
the same equations by an obvious replacement: m —>■ rrii, e — > ^j, x — > Xj, z = g, s . 

However, on this way there still exists a problem for the full baryon wavefunction. 
We did not exclude clearly the baryon center-of-mass motion from the wavefunction and 
provided only the zero average value of the baryon 3- momentum, see fll2j) . This means 
that the center-of-mass oscillates near zero point and bears a nonzero kinetic energy, 
which must be excluded from the baryon physical mass. The simplest way to do this is 
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to define the baryon mass in accordance with the Klein-Gordon equation: 

= E^- (p2), (25) 

where P = ^ pi is the total momentum of the three-quark system, and E = Y^Si is the 
total energy. 

As we consider baryons in the leading valence approximation with quarks in the 5"- 
wave orbitals, it is obvious that 



(PiPi) = for i 7^ j, and (p^) = {Xi\p^\Xi) 



(26) 



where Xi the eigenfunction of the system of equations (PTjl. Thus 

Hxi = (ap + (3-mi + i- (^-^^ + p ■ ar^ = e^i (27) 
and the averaged momentum squared can be expressed via two simple integrals: 



{Xi\P \Xi) = {Xi\ocp\oLp\xi) = ( Xi 



Xi 



r 



Ei- (3 ■ nii-l 

Xi)-'^lx 



(3 ■ ar 



Xi 



C 



+ - I (mj + ar) ■ (3 



Xi 



j dr {G\r) + F\r)) (^e, + ^) + + j - 

-2 J dr {G\r) - F\r)) (^e^ + (m^ + ar) (28) 



Now the baryon mass can be calculated as follows: 



M 



(29) 



As before, we neglect here the light quark mass, so, to make the numerical estimations, 
one needs to fix two dimensionful model parameters: the string tension a and the strange 
quark mass m^. Being guided by our previous calculations of baryon spectrum [Hj, we 
use the following values: 



a = 0.18 GeV^ = 210 MeV 



(30) 



In the Table E] we present the values of masses of ground-state baryons, calculated 
in this approach in comparison with experimental ones. One can see a fine agreement 
between the data. 

This work was supported in part by the Grant for support of Leading Scientific Schools 
# 843.2006.2, by the State Contract # 02.445.11.7424, and in part by the RFBR grants 
## 05-02-17869, 06-02-17120, 06-02-17012. 
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Table 2: Ground-state baryon masses, calculated in the Dirac orbital model, in compar- 
ison with the experimental values (in MeV) 



Baryon 


A 


E* 




n 


Theor. 


1233 


1381 


1527 


1672 


Exp. 


1232 


1383 


1532 


1672 
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